Decoherence-induced geometric phase in a multilevel atomic system 
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We consider the STIRAP process in a three-level atom. Viewed as a closed system, no geometric 
phase is acquired. But in the presence of spontaneous emission and/or collisional relaxation we show 
numerically that a non- vanishing, purely real, geometric phase is acquired during STIRAP, whose 
magnitude grows with the decay rates. Rather than viewing this decoherence-induced geometric 
phase as a nuisance, it can be considered an example of "beneficial decoherence" : the environment 
provides a mechanism for the generation of geometric phases which would otherwise require an extra 
experimental control knob. 



I. INTRODUCTION 

Berry observed that quantum systems may retain a 
memory of their motion in Hilbert space through the 
acquisition of geometric phases Remarkably, these 
phase factors depend only on the geometry of the path 
traversed by the system during its evolution. Soon af- 
ter this discovery, geometric phases became a subject of 
intense theoretical and experimental studies 0. In re- 
cent years, renewed interest has arisen in the study of 
geometric phases in connection with quantum informa- 
tion processing [H, 0]. Indeed, geometric, or holonomic 
quantum computation (QC) may be useful in achiev- 
ing fault tolerance, since the geometric character of the 
phase provides protection against certain classes of er- 
rors || 0, 0, [|| • However, a comprehensive investigation 
in this direction requires a generalization of the concept 
of geometric phases to the domain of open quantum sys- 
tems, i.e., quantum systems which may decohere due to 
their interaction with an external environment. 

Here we consider the following basic question: 

Is it possible for the environment to induce 
a geometric phase where there is none if the 
system is treated as closed? 

Apart from its fundamental nature, this question is of 
obvious practical importance to holonomic QC, since if 
the anwer is affirmative the corresponding open-system 
geometric phase can either be detrimental (if it causes a 
deviation from the intended value) or beneficial, in the 
sense that the environment is acting as an amplifier for, 
or even generator of, the geometric phase. 

Geometric phases in open systems, and more recently 
their applications in holonomic QC, have been consid- 
ered in a number of works, since the late 1980's. The 
first, phenomenological approach to the subject used 
the Schrodinger equation with non-Hermitian Hamilto- 
nians 0, [l(| • While a consistent non-Hermitian Hamil- 
tonian description of an open system in general requires 
the theory of stochastic Schrodinger equations [TTj] , this 
phenomenological approach for the first time indicated 
that complex Abelian geometric phases should appear 



for systems undergoing cyclic evolution. In Refs. [12|, [lj, 
HH EH EH E3 > geometric phases acquired by the density 
operator were analyzed for various explicit models within 
a master equation approach. In Refs. [8l. Il8|. the quan- 
tum jumps method was employed to provide a definition 
of geometric phases in Markovian open systems (related 
difficulties with stochastic unravellings have been pointed 
out in Ref. In another approach the density oper- 

ator, expressed in its eigenbasis, was lifted to a purified 
state [20, [21J • In Ref. [22| , a formalism in terms of mean 
values of distributions was presented. An interferomet- 
ric approach for evaluating geometric phases for mixed 
states evolving unitarily was introduced in Ref. |23ll and 
extended to non-unitary evolution in Refs. 0, Ha|. This 
interferometric approach can also be considered from a 
purification point of view [23l . [25j . This multitude of 
different proposals revealed various interesting facets of 
the problem. Nevertheless, the concept of adiabatic ge- 
ometric phases in open systems remained unresolved in 
general, since most of these treatments did not employ 
an adiabatic approximation genuinely developed for open 
systems. Note that the applicability of the closed systems 
adiabatic approximation [261 ] to open systems problems 
is not a priori clear and should be justified on a case- by- 
case basis. Moreover, almost all of the previous works 
on open systems geometric phases were concerned with 
the Abelian (B erry phase) case. Exceptions are the very 
recent Refs. [M [23, [2^| , which discuss both non-adiabatic 
and adiabatic dynamics, but employ the standard adia- 
batic theorem for closed systems in the latter case. 

Recently, a fully self-consistent approach for both 
Abelian and non-Abelian adiabatic geometric phases in 
open systems was proposed by Sarandy and Lidar (SL) 
in Ref. . It applies to the very general class of sys- 
tems described by convolutionless master equations 13011 . 
SL made use of the formalism they developed in Ref. [31| 
for adiabaticity in open systems, which relies on the Jor- 
dan normal form of the relevant Liouville (or Lindblad) 
super-operator. The geometric phase was then defined 
in terms of the left and right eigenvectors of this super- 
operator. This definition is a natural generalization of the 
one given by Berry for a closed system, and was shown 
to have a proper closed system limit. The formalism was 
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FIG. 1: Three-level atomic configuration with degenerate 
ground state levels |1), |2), and excited state |3). The atom in- 
teracts with two resonant classical fields with time-dependent 
Rabi frequencies gi(t) (probe laser) and Q2 (t) (Stokes laser). 



illustrated in Ref. [29[ in the context of a spin system 
interacting with an adiabatically varying magnetic field. 

In order to address the basic question posed above, 
we study here the adiabatic geometric phase in a multi- 
level atomic system using the SL formalism. Specifically, 
we consider the process of stimulated Raman adiabatic 
passage (STIRAP) [H, HH in a three-level atomic sys- 
tem in a A configuration. We analyze a version of STI- 
RAP where the closed system geometric phase is identi- 
cally zero. We then show that when spontaneous emis- 
sion and/or collisional relaxation are included, the same 
STIRAP process yields a non- vanishing geometric phase. 
This decoherence-induced geometric phase is an exam- 
ple of "beneficial decoherence" , where the environment 
performs a potentially useful task. This is conceptually 
similar to the phenomenon of decoherence-induced en- 
tanglement [34l. [35j. 

Since the SL formalism involves finding the Jordan nor- 
mal form of a general matrix, which is an analytically 
difficult problem, we developed a numerically stable pro- 
gram to find the Jordan form of any complex square ma- 
trix and used it to find the geometric phase [361 ] . 

The structure of the paper is as follows. In Section HT1 
we briefly review the STIRAP process in a closed three 
level system in the A configuration, and the correspond- 
ing calculation of the (vanishing) geometric phase. In 
Section IIIII we revisit this problem in the open system 
setting and derive the solution of the STIRAP model. 
Our numerical results, along with a detailed analysis of 
the geometric phase, are presented in Section [TV] We 
conclude in Section W\ 



II. GEOMETRIC PHASE UNDER STIRAP: THE 
CLOSED SYSTEM CASE 

We consider the process of stimulated Raman adia- 
batic passage (STIRAP) [33] in a three-level system in 
the A configuration, as shown in Fig. [T] In this process, 
the initial atomic population in level |1) is completely 
transferred to level |2), while the pulses are applied in a 
"counterintuitive" sequence. The intermediate level |3) 
does not become substantially populated. The interac- 



tion picture Hamiltonian in one-photon resonance can be 
written in the rotating- wave approximation as follows: 



H = gi (t)\3)(l\+g 2 (t)\3)(2\+h., 



(1) 



where the real functions gi are the time-dependent Rabi 
frequencies of the two laser pulses, interacting respec- 
tively with the transitions \i) *-> |3) (i € 1, 2). The eigen- 
values of H are given by 



E = 0, E ± = ± 



9l 
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(2) 



and the respective eigenvectors are given by 



|0) = cos(0)|l) -sin(0)|2), 

|+) = sin(0)sin(0)|l) + cos(0) sin(0)|2) +cos(0)|3), 
|-) = sin(0)cos(0)|l) + cos(0)cos(0)|2) - siru>)|3), 

(3) 

where tan(0) = 31/(72- Thus the time-dependence of the 
eigenf unctions is parameterized by that of 0. In principle 
the gi's can be complex valued, which gives rise to a 
controllable phase <p Q. Here we work with real valued 
gi& and set <f> — 7r/4 for the remainder of this work. The 
state |0) is a dark state, i.e., it has eigenvalue 0. 

We choose a Gaussian time-dependent profile for the 
control pulses: 



0i(f) = 5 oie- (t - to)2/r2 , g 2 (t) = g 02 e 



(4) 



where goi and 1702 are the pulse amplitudes, and to is the 
time-delay between the pulses, with pulse g 2 preceding 
pulse <?i . All time-scales are normalized in terms of the 
pulse-width r. The closed-system adiabaticity condition 
is satisfied provided to ~ r and [HJ 



dt 



< 1 Vt. 



(•5) 



In this limit, the evolution of the system strictly follows 
the evolution of either of the adiabatic states. Due to 
the ordering of the pulses as in (|4]), the atom initially in 
the level |1) is prepared in the adiabatic state |0). The 
population in level |1) is then completely transferred to 
level 1 2) adiabatically, following the evolution of state |0) 
under the action of the pulses (|4]) . Note that as the sys- 
tem follows the evolution of the state |0) in the adiabatic 
limit and the excited level |3) does not contribute to |0), 
the traditional view of the process is that it remains un- 
affected by spontaneous emission. Below we will show 
how this view must be modified in a consistent treat- 
ment of the process as evolution of an open system. In 
addition, incoherent processes such as dephasing of the 
ground state levels will affect the population transfer pro- 
cess. 

The geometric phases acquired by each adiabatic state 
|n), as acquired during the evolution between t and t 
can be easily calculated from [l[ 



n = i jf dt'{n\ — \n) 



(6) 
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consider the adiabatic eigenstates of Eq. ([3]). Then 
= |(sinW(l|+cosW<2|+cos(0)<3|)(cosW|l)- 
sm(0)|2)) - 0, and (-|$R = |(sm(0)<l| + cos(0)<2| - 
sin(0)(3|)(cos(6»)|l) - sin(0)|2)) = 0, irrespective of 
the value of <j>. Also, (0|^|0) = (cos(0)(l| - 
sin(0)(2|)(-sin(0)|l) - cos(0)|2)) = 0. Thus the STI- 
RAP process under consideration does not give rise to a 
closed-system geometric phase. We note that the analysis 
above is a special case of the four-level model considered 
in Ref. (H- 



III. GEOMETRIC PHASE UNDER STIRAP: 
THE OPEN SYSTEM CASE 



FIG. 2: A closed curve in the (51,52) parameter space, for 
to — t = goi = go2 = 1. At t = —00 the curve is at the 
origin, then rises steeply, and eventually returns to the origin 
at t — +00. 



In terms of a vector R(t) in parameter space undergoing 
cyclic evolution, this phase can be rewritten as 



f3 n = i <b( 



[n (R)\-^\n(R)) ■ dR. 
dR 



(7) 



In the case of the three-level system depicted in Fig. [T] 
the parameter space is defined by 31 (t) and 32 (t), i.e., 



fin = i^2 j>{n{9x,92)\-Q-\n{gi,g2))dgj. 

■j = l y 3 



(8) 



We consider a cyclic evolution in this parameter space, 
which takes place as t varies from —00 to +00, i.e., 
gi (-00) = 52 (-00) = 0, 31(00) = g 2 (oo) = 0. This 
is shown in Fig. [5J One can also parametrize the time- 
dependence of the pulses in terms of an angle 0, using 
(01 , such that 



tan0(i) 



SiW = £oi e ( 2ttQ - t g)/ T 

92{t) 302 



(9) 



Then as t varies from —00 to +00 we have that tan0(t) 
varies from to 00, and hence 8(t) varies from to 7r/2. 
Changing variables in Eq. (JT]), the geometric phase be- 
comes in our case: 



(3n 



* \n{d)\±\n{e))de. 



(10) 



Note that the relevant parameter space for our problem 
is that with coordinates (31,32), not (0, ip) of Eq. ([3]); 
indeed, Eq. (fTO]) does not even describe a cycle in the 
(8, ip) space, whereas the expression §8$ along with Fig. [2] 
show clearly that there is a cycle in the (31,32) space. 

Let us now show that the geometric phase van- 
ishes for all three adiabatic eigenstates \n) of ([3]), 
because the integrand (n(9)\jg\n(9)) = 0. Indeed, 



A. The model 

We now analyze the effect on the geometric phase of in- 
teraction of the atomic system with a bath causing spon- 
taneous emission and collisional relaxation. We describe 
these processes in the Markovian limit for the bath, us- 
ing time-independent Lindblad operators and neglecting 
Lamb and Stark shift contributions [3(3]. Thus the time- 
dependence appears only in the control Hamiltonian H 
[Eq. ([I])]. and the evolution of the system density matrix 
p is given by the Lindblad equation (in ft = 1 units): 

dp/dt = Lp = —i[H, p] + Cp, 
1 " 

c P = -^(2r l/9 rJ- P r|r l -rlr 4 p), (11) 

1=1 

where the dissipator C describes the incoherent processes, 
arising from system-bath interaction. We include spon- 
taneous emission from level |3) at rates 713 and 723 via 
Lindblad operators 



Ti =7i 3 |l)(3| ( r 2 = 7 23 |2)(3|. 



(12) 



We also include collisional relaxation between levels |1) 
and 1 2) at rates 712 and 721 via Lindblad operators 

r 3 = 7 i2|l)(2| ,r 4 = 72i|2)(l|. (13) 



B. Review of open systems geometric phase 

To see how a geometric phase can be associated with 
the master equation evolution, we follow Ref. [29[ and 
write the master equation as 



dp/dt = L[R(t)]p(t), 



(14) 



where L depends on time only through a set of parame- 
ters R(t) = R. These parameters will undergo adiabatic 
cyclic evolution in our problem. 

In the superoperator formalism, the density matrix 
for a quantum state in a Z?-dimensional Hilbert space 
is represented by a D 2 -dimensional "coherence vector" 
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\p)) = (Pii P2, ■ • • , Pd 2 Y (where t denotes the trans- 
pose) and the Lindblad superoperator L becomes a 
D 2 x D 2 -dimensional supermatrix [33], so that the mas- 
ter equation Q14p can be written as linear vector equa- 
tion in _D 2 -dimensional Hilbert-Schmidt space, in the 
form d\p))/dt = L[R(t)]\p)) . Such a representation can 
be generated, e.g., by introducing a basis of Hermitian, 
trace-orthogonal, and traceless operators [e.g., the D- 
dimensional irreducible representation of the generators 
of su(-D)], whence the pi are the expansion coefficients 
of p in this basis (37[, with pi the coefficient of I (the 
identity matrix). 

The master equation generates a non-unitary evolu- 
tion since L is non-Hermitian. In fact, L need not even 
be a normal operator (jJL ^ LlJ). Therefore L is gen- 
erally not diagonalizable, i.e., it does not possess a com- 
plete set of linearly independent eigenvectors. Equiva- 
lently, it cannot be put into diagonal form via a simi- 
larity transformation. However, one can always apply 
a similarity transformation S to L which puts it into 
the (block-diagonal) Jordan canonical form [38[ , namely, 
Lj = S~ 1 LS. The Jordan form Lj of a D 2 x D 2 matrix 
L is a direct sum of blocks of the form Lj = Ja (a 
enumerates Jordan blocks), where m < D 2 is the number 
of linearly independent eigenvectors of L, n a = D 2 

where n a = dim J a is the dimension of the ath Jor- 
dan block, and J a = X a I Ua + K a where X a is the ath 
(generally complex-valued) Lindblad-Jordan (LJ) eigen- 
value of L (obtained as roots of the characteristic polyno- 
mial), I na is the n a x n a dimensional identity matrix, and 
K a is a nilpotent matrix with elements (K a )ij — 
(l's above the main diagonal), where 6 is the Kro- 
necker symbol. Since the sets of left and right eigen- 
vectors of L are incomplete (they do not span the vec- 
tor space), they must be completed to form a basis. In- 
stantaneous right {\V { ^[R(t)]))} and left {((^ [R(t)}\} 
bi-orthonormal bases in Hilbert-Schmidt space can al- 
ways be systematically constructed by adding n a — 1 
new orthonormal vectors to the ath left or right eigen- 
vector, such that they obey the orthonormality condi- 
tion {{£a = SapS 1 ! [3l|. Here superscripts enu- 
merate basis states inside a given Jordan block (i,j S 
{0,...,n a -l}). When L is diagonalizable, {\V { ^ [R(t)]))} 

and {{{£& [R(t)]\} are simply the bases of right and left 
eigenvectors of L, respectively. If L is not diagonaliz- 
able, these right and left bases can be constructed by 
suitably completing the set of right and left eigenvec- 
tors of L (which can be identified with columns of S and 
S T , respectively, associated with distinct eigenvalues A Q ). 
Then for all times t 

L|X#>)) = p(J-V))+\ a \vW)), 
«40|L = ((£^\ + X a ((£^\, 

so that the {\V { i ] ))} and {((^1} preserve the Jordan 
block structure (see Appendix A of Ref. [2^| for a de- 
tailed discussion of these issues). 



In order to define geometric phases in open systems, 
the coherence vector is expanded in the instantaneous 
right vector basis {\V ( ^ [R(t)]))} as 

m rip— 1 

|P(*)» = E E pf(t)eti X ^ dt ' \vf[R(t)])), (15) 

/3=1 j=0 

where the dynamical phase exp^' \p(t')dt'\ is explicitly 
factored out. The coefficients {p^ (t)} play the role of 
"geometric" (non-dynamical) amplitudes. We assume 
that the open system is in the adiabatic regime, i.e., Jor- 
dan blocks associated to distinct eigenvalues evolve in a 
decoupled manner 31]. Then: 

np—l 
0\\ =\ a j=0 

Note that, due to the restriction Xp — X a , the dynamical 
phase has disappeared. 

A condition on the total evolution time, which allows 
for the neglect of coupling between Jordan blocks used in 
deriving Eq. (| 16[) . was given in Ref. [3lj ] . This condition 
generalizes the standard closed-system adiabaticity con- 
dition [26(, from which Eq. ([5]) is derived. Nevertheless, 
we have used the simpler condition ([5]) in our simulations 
below, as it is rather accurate in the present open system 
case. 

For closed systems, Abelian geometric phases are as- 
sociated with non-degenerate levels of the Hamiltonian, 
while non-Abelian phases appear in the case of degen- 
eracy. In the latter subspace of the Hilbert 
space acquires a geometric phase which is given by a 
matrix rather than a scalar. For open systems, one- 
dimensional Jordan blocks are associated with Abelian 
geometric phases in the absence of degeneracy, or with 
non-Abelian geometric phases in case of degeneracy. 
Multi-dimensional Jordan blocks are always tied to a 
non-Abelian phase (29j . 

1. The Abelian case: generalized Berry phase 

Consider the simple case of a non-degenerate one- 
dimensional Jordan block (a block that that is a 1 x 1 
submatrix containing an eigenvalue of L). In this case, 
the absence of degeneracy implies in Eq. (fit?)) that = 
X a a = (3 (non-degenerate blocks). Moreover, since 
the blocks are assumed to be one-dimensional we have 
n a = 1, which allows for removal of the upper indices 
in Eq. (|16p . resulting in p a = — p a ((£ a \'Da))- The solu- 
tion of this equation is p a {t) — p a (0) exp[i/3 a {t)], with 
I3 a (t) = i $l((£ a (t')\V a (t')))dt' . For a cyclic evolution in 
parameter space along a closed curve C, one then obtains 
the Abelian geometric phase associated with the Jordan 
block a 0: 

p a {C) =i £((£ a (R)\^\V a (R))) ■ dR. (17) 
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This expression for the geometric phase bears clear sim- 
ilarity to the original Berry formula, Eq. (0. Note that 
in general (3 a (C) can be complex, since ((£ a \ and \V a )) 
are not related by transpose conjugation. Thus, the geo- 
metric phase may have real and imaginary contributions, 
the latter affecting the visibility of the phase. As shown 
in Ref. [2!|, the expression above for (3 a {C) satisfies a 
number of desirable properties: it is geometric (i.e., de- 
pends only on the path traversed in parameter space), it 
is gauge invariant (i.e., one cannot modify the geometric 
phase by redefining ((£ a \ or \T) a )) via multiplication of 
one of them by a complex factor); it has the proper closed 
system limit (if the interaction with the bath vanishes, 
(3 a (C) reduces to the usual difference of geometric phases 
acquired by the density operator in the closed case). 



2. The non-Abelian case 

Ref. [2!| also derived the non-Abelian open systems ge- 
ometric phase, for the case of degenerate one-dimensional 



Jordan blocks. A non-Abelian geometric phase in fact 
arises in our STIRAP model when the spontaneous emis- 
sion rates are equal. However, we shall not treat this case 
in the present paper. 



C. Solution of the STIRAP model 

Returning to the STIRAP model, let us represent the 
density matrix p in terms of the coherence vector v as 



1 

N 



N(N - 1) 



(18) 



where the f2 Q are the Gell-Mann matrices [39(. Writ- 
ing the Lindblad equation p = Lp in the {tt a } basis, 
we obtain v = Lv, where v = |[1, s/3vi, ■ ■ ■ , v^us]* is 
a nine-component coherence vector. In the same basis 
we can express the Liouville operator L in the following 
form: 
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where we have used the Hamiltonian (fTJ). Here 

(713-723) > 



7+ 



(7i 2 3+7 2 2 3 )/2, 



7- 



7± 



(7 2 2±72i) I 
(20) 

and (71,2 are given in Eq. (H|). 

The eigenvalues and left and right eigenvectors of L 
can be found in terms of the parameters j±, 7j_ and 171,2, 
but the expressions are very complicated. The analytic 
determination of the corresponding left and right eigen- 
vectors is cumbersome, so instead we have used a numeri- 
cal procedure, which is based on the discussion presented 
in subsection IIII Bl above [36| . 

To exhibit some of the analytic structure, we temporar- 
ily make the further simplification that the spontaneous 
emission rates are equal: 713 = 723 = 7. This is the 
case, e.g., for D2 transitions in 23 Na [4fJ. In addition we 
assume temporarily that the collisional relaxation rates 
vanish: 712 = 721 = 0. With these simplifications L has 
the following three sets of eigenvalues (the ordering of 



subscripts is explained below): 



= {0,- 


-7 2 , 


2 Q 

7 3P 


Q 


!> 


iV3, Q 
2 ^3P 


6P 




Q 


!> 


*V3, Q 
2 ^3P 


6P 





A 9 = (- 7 2 + 

A 2 = A 3 = i(- 7 2 + iVQY, A 7 = A 8 = i(- 7 2 - iVQ) 

(21) 



where 



P = (x + v^ 2 + (Q/3) 3 
Q = 4(. 9 2 +. 9 2 )- 7 4 , 



1/3 



x = 7 2 (.9i+.9 2 2 ) 



(22) 



and the last set of four eigenvalues appears in two de- 
generate pairs. Because of this, the corresponding open 
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systems geometric phase is non-Abelian (recall the dis- 
cussion above), but we do not consider this case here. 

In the closed system limit (7 — > 0) L becomes — i[H, ■} 
and its eigenvalues are e nm — i{E n — E m ) (n,m G 
0, +,— ), where E n , m are the eigenvalues of the control 
Hamiltonian H as given in Eq. ([2]). The grouping in 
Eq. (j2"Tj) represents this limit in the following sense: 

A4, A5, A@ — > e nn = 0; 

Ai -> e + _ = 2i^g\+gl; 

A9 e_ + = -2iJgf+g$; 



(A 2 -► e+o) = (A 3 -> e _) = i\/g{ + g\\ 



(A 7 -» e 0+ ) = (A 8 -» e_ ) - -i\}g\ + gl (23) 

The corresponding eigenvectors of L reduce to |n)(m|. 
The subscripts of the A a represents the ordering of the 
eigenvalues in the closed system limit. We find that the 
degeneracy leading to a non-Abelian geometric open sys- 
tem phase appears only when 713 = 723 and 712 = 721 = 
0, or in the closed system limit. 

By a coordinate transformation from the control fields 
<7i ; 2 to the angle 9 = arctan(<7i/<72) we have, similarly to 
the closed system case, from the generalized geometric 
phase formula Eq. (fT7|) : 



a = i 



-it/2 



d 0((£ a \-\V Q )) 



(24) 



This expression for the phase associated with the ath 
eigenvector of L yields, in the closed system limit, not 
the absolute phase of each of the adiabatic eigenstates 
of the system Hamiltonian, but rather their phase differ- 
ences. This is natural as only a phase difference is an 
experimentally measurable quantity. 
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FIG. 3: Spontaneous emission without collisional relaxation: 
Variation of Re(/3i) (thin lines) and — Re(/3g) (thick lines), in 
units of 2-7T, with respect to 723T, for three different ratios 
between 723 and 713. The other parameters used are: 712 = 
721 = 0, t = 4r/3, g 01 r = g 02 r = 15, and r = 1. 




FIG. 4: Effect of reversing the order of the control pulses. 
Now to = — 4t/3 and all other parameters are the same as in 
Fig- El 



IV. RESULTS AND DISCUSSION 

We plot the real part of the open system Abelian ge- 
ometric phase, i.e., Eq. |24]) . for various combinations of 
the spontaneous emission and collisional relaxation rates 
in Figs. [3 1101 The main finding is that the answer to the 
question we posed in the introduction, "Is it possible for 
the environment to induce a geometric phase where there 
is none if the system is treated as closed?", is affirmative. 
Indeed, a glance at Figs. [3TTU1 reveals that the geometric 
phase is non-zero, and in fact increases with the decay 
rates. Moreover, we find that the imaginary part of the 
geometric phase is always zero to within our numerical 
accuracy, implying that the visibility of the geometric 
phase is unaffected in the present case by the interaction 
with the environment. 

In Fig. [3] we show the real part of the open system 
geometric phases [3\ and — /3g, for the case when colli- 
sional relaxation vanishes (712 = 721 = 0) and there is 




Y12- 

FIG. 5: Collisional relaxation without spontaneous emission: 
Variation of Re(/3i) (thin lines) and -Re(/3g) (thick lines) in 
units of 2tt with respect to 712T for different ratios between 
712 and 721. The other parameters used are: 713 = 723 = 0, 
to — 4r/3, goiT = <?02T = 15, and t = 1. 
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FIG. 6: Spontaneous emission without collisional relaxation: 
Variation of real part of phases in units of 2ir with respect to 
723T for 723 = 2713. Other parameters: to = 4r/3, goir = 
go2T — 15, and r = 1. 
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FIG. 8: Collisional relaxation without spontaneous emission. 
Variation of real part of phases in units of 2ir with respect to 
712T for 723 = 713 = and 712 = 2721. Other parameters as 
in Fig. © 




FIG. 7: Spontaneous emission without collisional relaxation. FIG. 9: Collisional relaxation without spontaneous emission. 
Parameters the same as in Fig. [5] except that 723 = 713/2. Parameters the same as in Fig. |S] except that 712 = 721/2. 



only spontaneous emission (713,723 7^ 0). Clearly, the 
phases increase monotonically with the emission rates. 
It is interesting to note that (to within our numerical ac- 
curacy) Re/?i = — Re/39 when 713 = 723. Recalling that 
Ai -> i(E + - EJ) and A 9 -> i(E- - E+) [Eq. (|25))]. this 
symmetry can be traced back to the difference between 
the adiabatic eigenstates |+) and |— ), which differ only in 
the sign of the coefficient in front of the excited state |3) 
[recall Eq. © and that (j> = 7f/4]. When the spontaneous 
emission rates are equal this difference in sign between 
the (|3) component of the) states |+) and |— ) generates 
only a difference in sign between the corresponding geo- 
metric phases, but not in magnitude, i.e., f3\ = — /3g. 

We also note that, in spite of the symmetry between 
the states |1) and |2) in our model, there is an asym- 
metry between the curves 723 = 2713 and 723 — ^713 
in Fig. [3] for a given geometric phase, e.g., /3i. Indeed, 
one might have expected a symmetry under interchange 
of the indices 1 and 2, in the sense that, e.g., the points 
/M723 = 2) and /3i(723 = 1) on the curves 723 = 2713 and 
723 = 5713 respectively, should have overlapped. That 



this is not the case is because the order of the pulses 
172 (first) and g\ (second) breaks the symmetry between 
states |1) and |2). Indeed, Fig. 2] shows the results for 
/3i when the pulse order is reversed (now gi precedes (72), 
and as a consequence the order of the curves 723 = 2713 
and 723 = ^7i3 is now reversed as well. In other words, 
swapping the pulse order is equivalent to swapping the 
spontaneous emission rates 723 and 713. 

In Fig. [5] we show the real part of the open system 
geometric phases Pi and — /3g, for the case when spon- 
taneous emission vanishes (713 = 723 = 0) and there is 
only collisional relaxation (712,721 7^ 0). The results are 
qualitatively similar to those in Fig. [31 with the excep- 
tion that now Re/?i 7^ — Re/3g when 712 = 721- This 
symmetry breaking can be attributed to the fact that 
the collisional relaxation operators directly connect the 
states |1) and |2), whereas these states are only connected 
to second order under spontaneous emission and under 
the control Hamiltonian |1| . The other interesting differ- 
ence between Figs. [3] and [5] is that spontaneous emission 
only leads to larger values of the geometric phase than 
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FIG. 10: Spontaneous emission along with collisional relax- 
ation. Variation of the phases Re(/3i) (thin lines) and Re(f3g) 
in units of 2n with respect to 723T for 723 = e7i3 and three 
different ratios between 712 and 721. The other parameters 
used are: to = 4r/3, goir = g 2T = 15, and r = 1. 



FIG. 11: Natural log of the left-hand side of the adiabaticity 
condition © (i.e., r{t) = In (||£| /V#i(i) 2 + £ 2 (i) 2 )) for the 
parameters used in our simulations, i.e., to = 4r/3, <7oit = 
902T = 15, and r = 1. 



collisional relaxation only. 

In Figs. 05] we show the real part of the open system 
geometric phases 02,03 and —fa, (3%. All four phases 
involve the dark state |0) in the closed system limit. 
Recall from Eq. (|23p that the pair of eigenvalues A2.3 
becomes degenerate in the closed system limit, as does 
the pair Xj^. Figures [5] and [7] show the case of vanish- 
ing collisional relaxation but non-vanishing spontaneous 
emission. Figures [5] and [5] shows the opposite case of 
non-vanishing collisional relaxation but vanishing spon- 
taneous emission. It is interesting to observe that when 
there is only spontaneous emission, as in Figs. [6] and [7] 
the phases corresponding to degenerate eigenvalues in the 
closed system limit are identical to within our numerical 
precision. We do not have an intuitive explanation for 
this symmetry, which is absent when there is only colli- 
sional relaxation, as in Figs. [5] and [5] On the other hand, 
the asymmetry between Figs. |5] and [7] and between Figs. [5] 
and [9] can again be attributed to the symmetry breaking 
between levels |1) and |2), due to the time ordering of 
the control pulses. 

In Fig.[TO]we revisit (3\ and /3g, and turn on both spon- 
taneous emission and collisional relaxation, and consider 
irrational ratios of the various decay rates (in order to 
eliminate potential accidental degeneracies due to ratio- 
nal ratios). Indeed, all six curves are clearly separated, 
and judging by comparison to Fig. [3] the effect of in- 
cluding both decoherence mechanisms is to increase the 
magnitude of the geometric phases (i.e., the decoherence 
mechanisms cooperate rather than interfere). 

As a final note, we should point out that varying the 
Rabi frequencies g\ and 32 along a closed cycle in param- 
eter space (Fig. [2]), i.e., letting the time t vary from —00 
to +00 (which we implement in practice by integrating 
from — to 7r/2), is incompatible with Eq. ([5]) for all 
times t due to the finiteness of the parameters involved. 
Indeed, Fig. [TT] shows the left-hand side of Eq. ([5]) for 
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FIG. 12: The phase /3( is computed by integrating from 9i 
(corresponding to ti — — 3.06r) to 82 (corresponding to t^ = 
4.39t). In this time range the left-hand side of Eq. ([5} < 1/e. 
Plotted is the variation of the difference between f3\ and /3[ as 
a function of 723T for three different ratios between 723 and 
713. The other parameters are the same as in Fig. [3] 



the parameters we have used in our simulations. It is 
clear that the adiabaticity condition is satisfied only for 
— ti t ^ ^2, where t\ w —3 and £2 ~ 4.4. However, 
when we repeat our calculations of the open system ge- 
ometric phase with 6 varying between angles 6\ and 62 
corresponding to the times t\ and t% (i.e., not along a 
complete cycle in the <?i,<?2 parameter space), we find - 
as can be seen in Fig. [12]- that the effect on the geomet- 
ric phase is entirely negligible. This confirms that our 
choice of parameters satisfies the adiabatic limit for all 
practical purposes. 



9 



V. CONCLUSIONS 

Our study of STIRAP in an open three level quan- 
tum system reveals that the interaction with the environ- 
ment can endow a system with a geometric phase, where 
none existed without the interaction with the environ- 
ment. Mathematically, the vanishing geometric phase in 
the closed system case is attributable to the vanishing 
integrand in the Berry formula. In a certain sense this 
is easily understood as the result of having a geometric 
phase determined by only a single parameter (6) , whence 
no solid angle is traced out in parameter space. It would 
then be natural to conclude that, by including the in- 
teraction with the environment a non-zero solid angle 
is created, implying that in the presence of decoherence 
motion along an orthogonal direction in parameter space 
must have taken place. However, one must be careful in 
accepting this explanation, since in fact the polar angles 



6 and <f> do not properly describe the parameter space 
in our problem: indeed, 9 varies from to 7r/2 (while <f> 
is constant) and thus does not describe a closed path, 
while the correct parameter space is that defined by the 
pulse amplitudes gi and gi (see Fig. [2]). Thus a proper- 
explanation of the intriguing effect of an environmentally 
induced geometric phase is still lacking and will be un- 
dertaken in a future publication. Here we conjecture 
that this is due to the non-commutativity of the driv- 
ing Hamiltonian and the decohering processes we have 
considered. It should be possible to test this by using 
the quantum trajectories approach to the open systems 
geometric phase [8]. Another interesting open question 
is to what extent the finding presented here can be made 
useful in the context of holonomic quantum computing 
Q, i.e., whether can one constructively exploit the envi- 
ronmentally induced geometric phase for the generation 
of quantum logic gates. 
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